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Abstract. We describe an algorithm that, given a fc-tuple of permutations 
representing the monodromy of a rational map, constructs an arbitrarily pre- 
cise floating-point complex approximation of that map. 
J.^ We then explain how it has been used to study a problem in dynamical 

systems raised by Cui. 



1. Introduction 



Let S be a topological oriented 2-sphere. The branched coverings S — > S con- 
^ sidered in this article are all orientation preserving. Let Q := {Qi}i£i be a finite 

pH subset of S with / = Z/fcZ. In j9J, Hurwitz describes an elegant classification 

+^ of branched coverings S — > S with branch value set contained in Q in terms of 

^ admissible fc-tuples of permutations (<7j G &d)iei- A fc-tuple is admissible if: 

i " i • the permutations (<Ti)ig/ generate a transitive subgroup of &d, 

• o"i • <T2 • • • Cfc = id and 

• the cycle lengths satisfy the condition 

(1) Yl E (length(c) - 1) = 2d - 2. 

2£l c cycle 

lf~) Of <Ti 

See §l.l| for more details regarding the classification. 
CO It is easy, using a computer algebra system such as Gap [7j, to enumerate all 

admissible fc-tuples of permutations; it is an altogether different problem to con- 
struct an analytic model of a covering associated to a given admissible fc-tuple 
. . of permutations. The purpose of this note is to describe such an algorithm, its 

implementation, and an application to complex dynamics. 

1.1. Hurwitz's classification. Two branched coverings / :S->S and g : S — > S 
are equivalent if there is an orientation preserving homeomorphism h : S —¥ S such 
that g = f oh. Hurwitz's result is a classification of equivalence classes of coverings 
in this sense. 

Choose a basepoint * e S \ Q. For each i e /, choose a path ji joining * to Qi 
in S \ Q, in such a way that 

• the paths ji intersect only at *, 

• the paths (71, . ..,7k) are ordered cyclically counterclockwise around *. 
The fundamental group G = tt 1 (S \ Q,*) is generated by paths ji that follow 
7i, wind once counterclockwise around Qi, and return to * along 7^. It has the 
presentation 

G = (% i G 1 1 71 • 72 • • • % = id). 
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Let / : (S, C) — > (S, Q) be a covering branched over Q. Number . . . , *d} 
the /-preimages of *. Then, for each i E I and each m G {1, ...,d}, the path 7, 
lifts to a path starting at * m and ending at *„ for some n —: <Ji{m). This defines a 
permutation &i for each i £ I. Note that the A:-tuple (o~i)i e i is admissible: 

• since S \ C is connected, the group ((7^) is transitive on {1, . . . , d}; 

• since 71 • 72 • • • 7fc = id, we have that <t\ ■ 02 • • • = id; 

• computing the Euler characteristic of S \ C via the Riemann-Hurwitz for- 
mula yields ([I]) . 

Conversely, let {ai)i^i be an admissible /c-tuple of permutations. Define a 
branched covering as follows: start with d disjoint copies of S, cut open along 
the paths 7^. If o~i{m) — n, glue the right boundary of 7, on m-th sphere to the left 
boundary of 7* on the n-th sphere. This defines a covering with branch value set 
contained in Q. It is connected because (a*) is transitive on {1, . . . ,d}. The Euler 
characteristic of the cover is 2, because of ([I]) and the Riemann-Hurwitz formula; 
so it is a sphere. 

The fc-tuple (o~i)i£i must be considered up to diagonal conjugation by &d, which 
amounts to numbering the spheres differently. The constructions above then define 
a bijection between equivalence classes of branched coverings and equivalence classes 
of appropriate /c-tuples of permutations. 

A coarser equivalence relation on coverings has also been considered, but is not 
the main focus of this article: two coverings /, g : S — > S are Hurwitz equivalent if 
there exist homeomorphisms ho, h\ : S — > S with / o hi = ho ° <?• Hurwitz classes of 
coverings may also be classified by fc-tuples of permutations; namely, by the orbits 
on appropriate fc-tuples of the symmetric group &d (acting as above) and the pure 
braid group on k strings. The latter group's generators act by conjugating, for any 
two consecutive points Qi, Qi+i in Q, the permutations <7; and 0^+ 1 by OiOi+x. This 
amounts to changing the "spider" Ji by twisting the legs 7i and 7i+i around 
each other. 

1.2. Analytic models. Assume now Q C P X (C) and that / : S \ C ->• P X (C) \ Q 
is a covering map. Then, / defines holomorphic charts on S \ C and it is not 
difficult to see that the points in C are removable singularities: we denote by S/ 
the corresponding Riemann surface. By the Uniformization Theorem, there is a 
conformal homeomorphism 0/ : S/ — > P 1 (C). The map F := / o (fij 1 : P X (C) — > 
P 1 (C) is a holomorphic branched covering, i.e., a rational map. Assume g = / o 
h : S P X (C) for some homeomorphism h : S S. Let 4> g : S 9 ->• P X (C) 
be a conformal homeomorphism and set G := g o tp^ 1 : P 1 (C) — > P 1 (C) be the 
corresponding rational map. Then, H = <f)f o (f)^ 1 : P 1 (C) — > P 1 (C) is a Mobius 
transformation and F = G o H . 

Therefore, up to precomposition by a Mobius transformation, the rational map 
F only depends on the equivalence class of covering / : (S,C) — !• (P 1 (C),Q). We 
say that / is an analytic model. 

1.3. Dynamics. If = 3, then we may assume Qi = 00, Q2 = and Q3 = 1 
within P 1 (C). Furthermore, precomposing / by an appropriate Mobius transforma- 
tion, we may also assume that Q c Then / : P 1 (C) — > P 1 (C) becomes a 
post-critically finite map: a self-map of the sphere all of whose critical points have 



finite orbits. Pilgrim linked in 13 the "dessin d'enfant" (the full preimage of the 



segment [0, 1]) of / with a dynamical invariant, its "Hubbard tree". 
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Post-critically finite maps may be considered as purely topological objects, namely 
branched self-coverings of the topological sphere S; a fundamental theorem of 
Thurston (see |6|) proves (except in few well- understood, low-complexity cases) 
that such a branch covering has at most one holomorphic realization, up to conju- 
gation by a Mobius transformation; furthermore, if #Q = 3 as above, then it has 
precisely one holomorphic realization. 

We describe in Section [8] a question by Cui in the theory of holomorphic dynam- 
ical systems, and give an explicit holomorphic realization of a topological map he 
constructed. 

This will also be our running example in the text. With Qi = oo, Q2 = and 
Q3 = 1, the permutations representing the map are 

a x = (1, 7, 11, 2)(3, 8)(4, 5)(6, 10) (9, 12, 13), 
(2) ( t 2 =(1,3,12,4)(5,9)(6,7)(10,13,11)(2,8), 
(73= (1, 5, 13, 6)(7, 10)(2, 3)(8, 11, 12) (4, 9). 

We seek a degree- 13 rational map /. 

In this specific example, the search can be made more feasible as follows. Setting 
all critical points as unknowns and eliminating is out of the question. With a little 
faith that the symmetry between 00, 0, 1 translates to /, let p{z) = 1/(1 — z) be 
the rotation permuting 00, 0, 1, and note that P 1 (C)/p is a sphere, branched at the 
two fixed points of p. If / descends to a map g on P 1 (C)/p, then (after change of 
variables) it has the form g{z) = z(p(z) / q(z)) 3 for degree-4 polynomials, such that 
g(z) = 1 + 0((z — l) 4 ) at z = 1, and such that 1 is the image of four other points 
with local degrees 3,2,2,2 respectively. We are grateful to Noam Elkies and Curt 
McMullen for having pointed out to us the feasibility of this approach. 

Nevertheless, we will show that our algorithm is strong enough to produce a 
solution even without exploiting the symmetry of the Hurwitz data. 

1.4. Simple cases. If = 2, then there is a unique solution represented, up to 
diagonal conjugation, by the pair of permutations 

a\ = (1,2, ...,d) and cr 2 = (d, . . . , 2, 1). 

If Qi = 00 and Qi = 0, an analytic model is f(z) = z d . 

However, the case #Q = 3 seems already as complicated as the general case, 
and has only been addressed in the literature for small d. Such maps are often 
called "dessins d'enfant", see (si; the corresponding combinatorial objects for the 
modular surface f)/ PSL 2 (Z) are called "Conway diagrams", see [l] §3.4]. Methods 
of constructing them are addressed, inter alia, in [3jj5]. 

In this section, we consider the case d < 3 which can completely be solved. If 
#Q = 2, then as we said above we may choose Q — {00, 0} and f(z) = z d . If 
d = 4^Q — 3 then we may choose Q = {00, 0,1}. Without loss of generality, we 
may assume that all points of Q are branched values, since otherwise we are reduced 
to the case #Q = 2. Up to permutation of the points in Q and the indices, the 
only possible triple of permutations is 

<j x = (1, 2, 3), ct 2 - (1, 2) and a 3 = (2, 3). 

To find an analytic model, we seek a rational map / of degree 3 such that 

00 h> 00, h4 and 1 h4 1. 
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This implies f(z) = 3z 2 — 2z 3 as the only realization. 

The next case we consider is d = 3 and #Q = 4. Using Mobius transformations, 
we may normalise Q to be {oo, 0, 1, w}. Up to conjugation in ©3 we may take the 
first permutations to be a\ = . . . = <7j = (1, 2). The condition that the permutations 
generate a transitive group imply that one of them is not (1, 2). Up to conjugation, 
we may assume that the first permutation which is not (1, 2) is c^+i = (2, 3). Since 
(Ji • (72 • C3 • (T4 = id, this gives four possibilities, namely, writing a = (a±, 02, 03, 04), 

g = ((1,2), (1,2), (2,3), (2,3)), 0= ((1,2), (2,3), (1,2), (1,3)), 

(7= ((1,2), (2,3), (1,3), (2,3)), 0= ((1,2), (2,3), (2,3), (1,2)). 

To find the corresponding /, assume without loss of generality that / maps 00 i— > 00, 
1 y 0, 1 <— > 1 and v 1— >■ w. This forces the map / to have the form 



2 a(*-l) + l 



/a(z) ~ Z (a + 2)(z-l) + l' 

for some parameter a subject to (a + l)(a — l) 3 = u>a(a + 2) 3 ; then (a + l)(a — 1) = 
va(a + 2). Since w 7^ 0, 1, the equation defining a in terms of ui has four distinct 
roots, leading to four candidate maps f a . There is a bijection between the maps / 
and the triples of permutations above, but no canonical one — it will depend on 
the specific choice of #Q generators % of 7Ti(P 1 (C) \ Q,*)- Note also that these 
four solutions are part of a single Hurwitz class. 

2. Overview of the algorithm 

We are given a list . . . ,cr& of permutations in ©^ with product <7i • • • U}. — 1, 
and points Qi,.. . ,Qk € P 1 (C). 

Let cti = (a^i, . . . , ai.£ { ) be the cycle lengths of er^; we have ctij = d for all 
i, and X)i ~ 1) = 2d — 2. In the first part of the algorithm, we enumerate 

all rational maps with critical values Qi, . . . , Qk such that the multiplicities of the 
preimages of Qi are an, . . . ,a^i v In the second part, we select the appropriate 
rational map among these candidates. 

The approach in the first part of the algorithm seems to originate in Malle ITT) ; 



see also 12 



For the sake of describing its workflow more clearly, the actual algorithm (de- 
scribed in the remainder of the text) has been slightly simplified. 

Normalization: Without loss of generality, we assume Q\ — 00, Q2 — and 
Qs = 1. We approximate the other Qi by Qi £ P 1 (Q). The rational map 
we seek will leave Qi, Q2, Q3 fixed. 

Finite field solution: We pick a prime p, such that the points Qi have re- 
alizations Qi £ P 1 (F p ). We then list all degree-d rational maps F over F p 
with poles and zeroes of multiplicities a± and oti respectively, and by brute 
force check for each F whether F — Q { has zeroes of multiplicities a, for all 
i > 3. Note that the rational map F is a solution to our original problem 
over F p . (If there are no solutions, we restart with a different prime p). 

p-adic solution: Write F = W2/W1 with W2 monic of degree d, and W\ 
of degree less than d. (In fact, we later write the denominator as \W\ 
with W\ monic. The present discussion uses a simplified notation.) For 
i > 3, let Wi = W2 — Q,W\ be the numerator of F — Q { . We compute 
high-precision p-adic approximations Qi of the Qi, and lift each Wj from 
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F p to a high-precision polynomial Wi over Z p , in such a manner that we 
have Wi = W2 — QiW\ + 0(p N ) for large N. This lifting can be done 
by Hensel's lemma, because by Corollary 3.3 the Jacobian of the system 
{Wi = W2 — QiWi} is invertible at a solution for almost every prime p. (If 
DF happens not to be invertible, we restart with a different prime). 
Algebraic solution: Using the lattice-reduction algorithm LLL [To], we find 
polynomials Wi over Q, with coefficients of small height (small degree and 
coefficients of minimal polynomial) that are close to Wi obtained at the 
previous step. Using exact arithmetic over Q, we check that the solution 
W2/W1 is correct. (If not, we either compute a finer p-adic approximation, 
or higher-degree algebraic number approximations, or we restart altogether 
with a larger prime). 
Complex solution: For each root a* € Q of Wj, given by its minimal 
polynomial over Q, we compute (to user-specified precision) all the com- 
plex roots Ci,j, s of its minimal polynomial. We still need to find out 
which of the Galois conjugates Ci,j, s belong to the same solution. For 
this, we also choose small, random integers m, n, compute the minimal 
polynomial of mcij + nci' t j', and find complex approximations di^^j'.t 
of its roots. We then pair together those roots {(k,j siCi'j'y) f° r which 
di,j,i',j',t ~ mcij.s + nci' _ji tS i for some t. By considering enough of these 
pairs we can stitch together a collection of compatible coordinate approxi- 
mations Cj lJlS embracing all 

We call C the collection of all roots Cij, s , and note that the rational map 
is determined by its zeroes and poles, and the normalization that Q\ = 1 
is fixed. Since {0,oo} C Q, these zeroes and poles are in C. 
The second step of the algorithm checks, by path lifting, that the monodromy 

around Qi is correct. For each of the Galois conjugate solutions (C, /) obtained in 

the first step, we do the following: 

Triangulate: We are given a floating-point approximation Q of Q. We com- 
pute a triangulation i? of P 1 (C) whose vertex set contains Q, and a trian- 
gulation c € of P 1 (C) whose vertex set contains C. For efficiency reasons, 
we use Delaunay triangulations, see ^7] We compute the dual triangulation 
JS^; it has one vertex per face of J2, and edges transverse to those of =2. 
We fix a vertex * G £2^ as our basepoint. 

Lift the triangulation: Let W denote the vertices of J2- 1 . For each w € W, 
we number arbitrarily w%, . . . , Wj, the /-preimages of w. 

For each edge e G ^ _L , going from w' to w" , we compute a permutation 
<T e G &d such that the /-lift of e starting at w[ ends at u/' There are 
two strategies for this, one is by subdividing appropriately the path e and 
playing "connect-the-dots" , the other uses more efficiently the triangulation 
<€. 

Read permutations: For each critical value Qi £ Q, let e(l), . . . , e(n) be 
the sequences of edges traversed by a path in J2- 1 that starts and ends in 
the basepoint *, and surrounds once counterclockwise the point Qi and no 
other vertex of Q. Compute the permutation a[ = ^ e {i) ' ' ' ^e(n)- 

Check: The data (C, /) are a valid solution to the Hurwitz problem if and 
only if there exists a permutation r € &d such that a[ = (<Ji) T . 
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2.1. Implementation. The fourth-named author has implemented the first part 
of the algorithm, mainly in C, and the second-named author has implemented the 
second part of the algorithm, mainly in GAP |7 . By far the most time-consuming 



part of the procedure is the search for a solution over a finite field. Example 4.9 



required approximately 15 Minutes on a desktop, 30-SPECint2006 computer. The 
code is maintained by the fourth-named author, and is available at 

https : / / github . com/ j akobkroeker/HMAC 



3. The space of rational maps 

We show, in this section, that (as soon as the prime p is sufficiently large) we 
may lift every F p -solution to Z p . This follows from the well known fact that the 
Hurwitz spaces are smooth. We could not find the precise statement we need in 
the literature, so we give a complete proof. 

Let d > 2 be an integer and denote by Ratd the space of rational maps of degree 
d, which may be identified with a Zariski open subset of P 2d+1 (C). 

Let k > 3 be an integer and let F : S —> S be a ramified covering branched 
over Q = {Qi, . . . ,Qk\- Note that, according to the Riemann-Hurwitz Formula, 
C := F^ 1 (Q) contains exactly (k — 2)d + 2 points. We write C = [j i Ci with 
Ci := F~ 1 {Q i } = {Ci t i, . . . , Cj^}, and for each j e {1, . . . , we let a^j be the 
local degree of F at Cij. 

Let be the smooth quasiprojective variety of injective maps q : Q —> P 1 (C). 
For q € 13, we use the notation % := q(Qi). Similarly, let £ be the smooth 
quasiprojective variety of injective maps c : C — >• P 1 (C). For c € £, we use the 
notation Cjj := c(Cjj). The quasiprojective variety 2) := £ x Q is smooth. We 
shall prove that the subvariety 

£ := {(c,q) G 2) | (3/ e Rat,,) (Vi,j) f(c hJ ) = Ql and deg c . ./ = a itj } 

is also smooth, and regularly parametrised: 

Proposition 3.1. The variety X is smooth of dimension k + 3, locally regularly 
parametrised by (q u . . . , q k , c^i, c 2 ,i, c 3) i). 

Observe that, for (c, q) £ X, there is a unique rational map / € Ratd such that 
f( c i,j) ~ Qi an d deg c . ./ = otij for all (i, j). Indeed, knowing a rational map above 
three points completely determines the rational map (it is even enough to know the 
full preimage of two points plus one preimage of a third point). 

Note that the group of Mobius transformations acts on £ and by postcompo- 
sition: 

(M,N) ■ (c,q) := (M o c,A^o q). 
The quotient space may be identified with 2)o : — x 0o with 
Co := {c G £ | c\ t \ — oo, C24 = and 03,1 = l} 

and 

0o := {q G I qi = 00, q 2 = and q 3 = l}. 
The projection 2J -> 2J/(PSL 2 (C) x PSL 2 (C)) ^ 2J is a submersion. 
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The action preserves X as indicated on the following commutative diagram: 
C< ^pi(C)— ^P^C) 



I 



JVo/oM" 



qc — !Upi( C )_?Upi(c). 

It is therefore enough to show that Xo := Ifl 2)o is a smooth subvariety of 2)o 
locally regularly parametrised by (<?4, • • • , 

We first write equations for Xo- To each (c, q, A) G £o x fl x C*, we associate a 
collection of monic polynomials (Wi)ie{i,--Mk} defined by 

WiO) = ]J( Z - ci,j) ai ' j and for i > 2 Wi(z) := ]J(z - c^) "'" 
i=2 j=i 

and a collection of rational maps (/i)ie{2,....fc} defined by 

t _ jgj 

Note that these are degree-d rational maps with poles of order ax,j at c\j. In 
addition, fi maps Ci.j to qi with local degree oeij. It follows that (c, q) G Xo if and 
only if there is a A € C* such that /j = /2 for alii G {3, . . . , A:}, that is, -F, = with 

(3) F i :=W i + \q i W 1 -W 2 . 

In that case, we use the notation 

f(c,q,\) '■= h = /3 = • • • = fk- 

In other words, consider the map 

J- := (F 3 , . . . , F k ) : 2) x C -> (C[z] dog<d ) fc - 2 . 

Then, (c, q) G X if and only if there is a A G C* such that .F(c, q, A) = 0. 

According to the following Lemma and the Implicit Function Theorem, the sub- 
variety of 2)o x C defined by the equation F = is smooth of dimension k — 3, 
locally regularly parametrised by ((74, . . . , q^)- It follows that its projection to the 
2)o component, namely Xo, is also smooth of dimension k — 3, locally regularly 
parametrised by (94, . . . , q k )- 

Lemma 3.2. If J-(c, q, A) = 0, then the derivative D^ q x^IF restricts to an isomor- 
phism T c € x {0} x T\C -> T (C[z] dcg<d ) fc - 2 . 



We postpone the proof of the lemma to Section |3.2| and mention immediately a 
corollary that we shall use later. For q G 0, let J-" q : £0 x C* — > (C[z]d e g<d) k ~ 2 be 
defined by 

J" q (c,A) := J"(c,q,A). 

Corollary 3.3. Assume that (c, q, A) is defined over Q with J-"(c, q, A) = 0. Then, 
for almost every prime p, the derivative DTn at (c, A) is invertible mod p. 

Proof. Since the point (c, q,A) is defined over Q, its coordinates may be written 
using integers, and reduced mod p. For all except finitely many values of p, the 
resulting reduction gives a genuine point, namely where the reductions of (c, q) are 
injective and the reduction of /( c ,q.A) has degree d. Since DT is invertible over Q, 
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t=0 



it may be written as a/N for an algebraic integer a and N € N; then the reduction 
modulo p of DJ- is invertible for all primes not dividing N. □ 

Before embarking in the proof of the Lemma, we first build up a description of 
the tangent space of Rat^. 

3.1. The tangent space to rational maps. Consider a rational map / e Rat,j. 

A tangent vector to / is 

1 ' dt 

for a holomorphic family of rational maps (ft)\t\<e with /o = /. 

For every z £ P 1 (C), the vector f(z) is a tangent vector in P 1 (C) at f(z); in 
other words, / is a section of the pullback bundle /*(TP 1 (C)). It can be pulled 
back to a vector field on P 1 (C), as 

V (z) = -(D z f)- 1 (f(z)), 

or, in coordinates, r](z) = —f(z)/f'(z). Therefore, r/(z) is a meromorphic vector 
field on P 1 (C), holomorphic away from critical points of /, and with a pole of order 
at most m at critical points of multiplicity m. 

Geometrically r/(z) is the movement at time t = of the point z t — ft (/(z)). 
This point ft can be followed away from critical points, by the Implicit Function 
Theorem. 

We consider now local perturbations of / at a critical point, i.e. we assume 
that the vector field / is given by a path ft = 4 > t°f° 4't~ 1 with <f>t,ipt analytic 
perturbations of the identity at the critical value and point v, c respectively of /. 

Let Ct denote the critical point of ft and let vt denote its critical value; then 
c t = ipt (c) and Vt = (fit (v) ■ Let c denote the motion vector of Ct , and let v denote 
the motion vector of Vt- Then c = i))(c) and v — (fi{v). Now / = <fi o f — Df o t/j, 
because (f>o — ipo = id- Therefore, 

n + ri> = i>- 

At v, the vector field <fi takes value v, the vector field r\ + f*<fi is holomorphic at c, 
and its constant term is c. If v = 0, then f*<fi is holomorphic near c and vanishes 
at c. 

Therefore, whenever we have a family of rational maps (ft) for which we can 
follow a critical point Ct and its associated critical value Vt with v = 0, the vector 
field rj is holomorphic near c and coincides with c at c. 



3.2. Proof of Lemma 3.2 For i € {3, . . . k}, we have 

Ft 



fi - f: 



2 



XW 1 

Recall that if T(c, q, A) = 0, then fi = fi for all i € {3, . . . , k}. We denote by / 
this common rational map of degree d. If in addition(c, q, A) belongs to the Kernel 
of DT at (c, q, A), then 

i _ i Fj Fj- (AWi + XWi) = 
/4 /2 AWi (AT^i) 2 

so /i = /2 f° r all i € {3, ... , fc}. We denote by / this common tangent vector to 
Ratd at / and by rj the corresponding meromorphic vector field on P 1 (C). 
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As (c,q,A) varies in 2)o x C*, the /j-preimages of the points q\ = oo and Qi 
vary holomorphically: they are the points C\j and Cij. According to the previous 
remark, we see that if q = then, for all i e {2, . . . , k}, the meromorphic vector 
held -q is holomorphic near c\j, coincides with t\j at c\j, and furthermore is 
holomorphic near Cij and coincides with c,j- at c,j. 

Thus, r) is a holomorphic vector field on the whole sphere P 1 (C) and coincides 
with dij at Cij. In particular, it vanishes at c^i = oo, c 2 .\ = and c 3il = 1. A 
holomorphic vector field with at least 3 zeroes globally vanishes. Therefore, n = 
and bij = for all In addition, for all i e {1, . . . , k}, we have that W% = 

and for all i G {3, . . . k}, we have that = F t = AtfaWi. This shows that A = 0. 

Let us summarize: if J"(c, q, A) =0 and if (c, 0, A) belongs to the Kernel of DF at 
(c, q, A), then c = and A = 0. So, the restriction of Du^^T to T c £ x {0} x T\C 
is injective. Since T c £o x {0} x T\C and T (C{z]d C g<d) k ~ 2 have the same dimension, 
that is (k — 2)d, this restriction is an isomorphism as required. 

4. Finding a solution in a finite field 

We describe in this section an efficient method of finding a rational function over 
a finite field with prescribed critical values and multiplicities. 

We start by recalling some facts about univariate polynomials over non-algebraically- 
closed fields k of arbitrary characteristic. For this we need some notation: 

Notation 4.1. An ordered sequence a — (a\, . . . , otk) with a\ > ■ ■ ■ > > and 
a i — d is called a partition of d. With the shorthand notation 

0" := (^_^£) 

/i times 

We can always write a — (a\, . . . , a^) = ((3^ 1 , . . . , /3% n ) with fix > ■ ■ ■ > f3 n and 
appropriate /i^. For example, 13 = 4 + 3 + 2 + 2 + 2 is written as a = (4, 3, 2, 2, 2) = 
(4\3\2 3 ). 

The partition a* defined by a* := #{«: «i > j} is called the dual partition of a. 
For example, a* = (5, 5, 2, 1). 

Let / £ k[x] be a degree-ii polynomial, let k = (x — Ci) € k[x] be its distinct 
linear factors over the algebraic closure of k, and let a* be their multiplicities, so 
that 

k 

f = U i r 

i=l 

Without restriction we can assume a.\ > ■ ■ ■ > and a = (a\,...,a.k) is a 
partition of d. In this situation we say that / is of shape a. 

If we write a — (a\, ... ,ctk) = {Pi 1 , ■ ■ ■ ,Pn n ) as above, we can write 

n 

f = Uft 

with dcg/j = fii and fi the product of those linear forms that have multiplicity 
In this situation the /j are coprimc. A 

Lemma 4.2. Let k be a field of characteristic p, let f <G k[x] be a univariate 
polynomial of shape a and write f = fliLi w ^ h e -V a « < P f or a ^ * 
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then 

k 

gcd(/ ) /o=ii z ? i " 1 - 

i=l 

Proof. We have 

/'=x>£Mnc-Ox>(n<iX 

i—l i—1 j^i 

This shows that ]X=i 1 divides the gcd. Assume now that there is another 
linear factor I in the gcd. Since the gcd divides / there exists an index s with 
I = l s . Since the gcd divides /' we have that I divides 

«=1 j^i 

Now all summands except for a^fr^.^ lj)l' a are divisible by I. Since a s and l' s are 

nonzero in k and h[x] respectively, it follows that / must divide Yljjti^j- This is 
impossible since the l s are pairwise coprime. □ 

Corollary 4.3. With the notations of the previous Lemma we have 

g cd(/, /',..., /( e >)= n ir~ e . 

i: O-i >e 

Proof. Lemma |4.2| and induction. □ 
Corollary 4.4. With the notations above let a* be the dual partition of a. Then 

a* e = deggcd(/, /',..., - deggcd(/, /',..., /«). 

Proof. We have 

^ g cd(/, /v.., = n jl ^> g - j jr^ n 
9e ' g cd(/, /',..., /^) ai> jv- e 

It follows that 

deg^e^deg J| Zj = #{i : on > e} = a*. □ 

i : ai>e 

Algorithm 4.5 (Compute the shape of a polynomial). 

Given: a polynomial / 6 

Return: the shape (ai, . . . , afe) of /. 

Write d = deg(/). For each e = 0, . . . , d, compute g e = gcd(/, /',..., /^)- For 
each e = 1, . . . , d define then a* = g e -\ — g e . Return the dual of the partition 
(a* l7 ...,a* d ). 

Proof of validity. This directly follows from Corollary |4.4| □ 



We may collect linear factors of the same multiplicity, so as to avoid field exten- 
sions: 

Corollary 4.6. With the notation of Lemma \4-.2\ choose j3\ > ■ ■ ■ > f3 n among the 
ai such that f = FJ™ =1 /f* with fj = FJ i . ai _g. h- Then the fi are defined over k. 
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Proof. Since the calculation of a gcd does not require field extensions, we have 

gcd(/,/',...,/( e ))ek[x]. 

By Corollary |4.3| we then have 



_ gcd(/, f fVW) II: ft> ft -l /f'~ ft+1 * 

»- gcd (/ , /',..., /<*>) = Uu = l\ fiek[x] 

so 

fj = e k[x] n k(x) = k[x}. □ 
9j+i 

We are now ready to describe our algorithm searching for rational maps over ¥ p . 

Algorithm 4.7 (Compute all rational maps over k with given shape above given 
points). 

Given: a finite field k, a list of points Q 1; . . . , Q k £ P 1 (k), an integer d, and a list of 
partitions (ai, . . . , a k ) of d with a,; = (a^i, . . . , on t e t ) satisfying J2i J2j( a i,j — 1) = 
2d -2 S 
Return: all rational maps over k of degree d such that every Q i has ii preimages 
with local degrees a^i, . . . , respectively. 

We choose a Mobius transformation M € PSL2(k) sending Q 1 to oo and Q 2 to 

0. 

We write each partition on in compacted form as on = (Z^' 1 , ■ • ■ , /Sf^" 4 )- 
We enumerate all £i-tuples of monic polynomials (fx,... , f ni ) with deg(/j) = 
fiij, and all ^-tuples of monic polynomials (51,.. ■ ,ffn 2 ) with deg(gj) = (Jqj. 
For each such pair of tuples, we compute 

Wi = l[ f- 1 '" and W 2 = JJ 
j'=l j'=l 



Using Algorithm 4.5 we filter those (WijW-j) such that the shape of Wi is ot\ 
and the shape of W% is «2 (this fails only if a pair fi,gj is not coprime). 

By computing their g.c.d., we filter those (Wi,!^) such that W\ and W% are 
coprime. 

For each i = 3, . . . , k, l et C k be the set of A € k such that the shape 
(computed using Algorithm 4.5) of Wi :— W2 — XM(Q i )Wi is an. We filter those 
rational maps for which |~L=3 ^ * s non-empty. 

We return all the rational maps M -1 o (W2/AW1), for all A € 0^=3^, that 
survived the filtering. 

Proof of validity. Let first / := A/ -1 o (W2/AW1) be a rational map returned by 
the algorithm. For i = 2, . . . , k, consider the rational map fi — M^ 1 o (W^/AWi). 
By the very definition of Wi (compare with (|3|), we have / = f 2 = • • • = fk- On the 
other hand, the /j-preimages of Q i are the zeroes of Wi, so they have multiplicities 

on. _ 

On the other hand, let / be a rational map such that every Q i has £i preimages 
with local degrees on t i, ■ ■ ■ , otij i respectively. Then, for every Mobius transforma- 
tion M sending Q 1 to 00 and Q 2 to 0, the rational map M o f will be of the form 
W2/XW1, for monic polynomials W\,W% of respective shapes «i,a 2 and a scalar 
A G k. Furthermore, by Corollary 4.6 both W% and W2 factor over k[x] into poly- 
nomials of degrees fiij and [i 2 j respectively. The rational map Mof — Q i will be of 
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the form Wi/XWi with Wi of shape cti\ therefore, that solution / will be returned 
by the algorithm. □ 



Remark 4.8. Algorithm 4.7 is the most computationally-intensive part of our 
procedure. Its performance is improved in the following ways: 

(1) If = 1 for some i,j, then we may assume, after permuting the shapes 
a;, that i = 1 so that W% contains a power of a linear factor f? %,i . Fixing 
the corresponding preimage of oo to be oo amounts to the choice /j = 1, so 
that the degree of W\ is actually d — fyj. This speeds up the search by a 
factor p. 

Similarly, if up to permutation of the indices there are more ji^j = 1, 
with i £ {1,2}, then the corresponding factors may be assumed to be x 
and x — 1. 

On the other hand, if all fiij > 2, then no normalization of the critical 
points may be assumed, and in particular oo should not be assumed to be 
a preimage of some 



(2) When using Corollary 4.4 one can detect a wrong shape already if gcd(/, /') 
or for that matter any gcd(/, . . . , f^ e ') with e < ol\ has the wrong degree. 
We stop the calculation of gcd's as soon as this happens. This speeds up 
the process by a factor of about ol\. Similarly, as soon as the intersection of 
the Aj already computed is empty, the pair (Wi,Wa) should be discarded. 

(3) If the largest is small enough (e.g. 7 or 8 in Fn) we can enumerate all 
monic irreducible homogeneous polynomials of degree < Hij and build the 
fi and gi out of them, while taking care that no irreducible piece is used 
twice. We can then omit checking shape and coprimeness of W% and W 2 as 
these conditions are then automatically satisfied. 

Example 4.9. Over Fn we searched for a rational map of shape (4,3,2,2,2), 
(4,3,2,2,2) and (4,3,2,2,2); we chose Q 1 — oo and Q 2 = 0, and didn't specify 
Q 3 , letting on the contrary the algorithm determine choose it for us. We found the 
solution 

W 2 x 4 {x + 3) 3 (x 3 -3s- 5) 2 
Wx ~ {x-5) 3 {x 3 + 3x 2 + 2x + 3) 2 ' 
It has indeed the desired shape as we have the following factorisation 

W 2 + 4Wi =(x- l) 4 {x - 3) 3 (x 3 - 2x - 3) 2 , 

which implies, for the choice Q 3 = 1, the value A = —4. 

5. Lifting a solution from F p to Z p 

The lift from F p to Z p is done using Hensel's lemma (namely, Newton's method 
in positive characteristic): 

Proposition 5.1 (Hensel's Lemma). Let T — {F\, . . . , F m ) be a vector of poly- 
nomials, with Fi e Z[a^i, . . . ,x m ], and let J = (g^-) be the Jacobian matrix of F. 
Assume that a = (a±, . . . , a rn ) G Z m satisfies 

F(a)=0 modp N , 

that J(a) is invertible modulo p N , and let J~ 1 (a) be an inverse modulo p N . Then 

F{a + bp N ) = modp 2N , 
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for 



2N 



Furthermore J(a + bp N ) is invertible modulo p 

Proof. F(a) is divisible by p N since F(a) = mod p N . Therefore b is well defined. 
We have 



F(a + p N b) = F(a) + bJ{a)p N mod p 
F(a) 



= F(a)- -^J-\a)J(a)p N mod p 2N 
= mod p 2N . 

The invertibility holds more generally. Let A and £ be matrices with AB = 1 
mod p N . We can then write 

AB = t+p N C modp 2N . 

In this situation we have 

A(B - p N BC) = 1 + p N C - p N ABC mod p 2N 
= 1 mod p 2N 

since ^4_B = 1 mod p N ; so B' = B ~ p N BC is an inverse to A modulo p 27V . □ 

Consider the following data: a ring k; a family of polynomials VFi € k[a;] of degree 
at most d, for i = l,...,fc, with factorisations W, = Yl^Li^ij 1 1 a parameter 
A e k x ; and a sequence of points Qi = oo, Q2, ■ ■ ■ ,Qk € P 1 (k). We say that they 
are coherent if Wi/XWi + Qi is independent of i = 2, . . . , k. 

We say that they are normalised if the following holds: the first three values 
Qi, Q2, Q3 are 00, 0, 1 respectively; and the first three preimages C\ t \, 62,1, 63,1 are 
also respectively 00, 0, 1. This means that we assume that W\ has degree d — ot\ t i, 
that W 2 ,i = x a2 '\ and that W 3<1 = (x - l)" 3 - 1 . 

Note that this assumption is not innocuous: it may well be that no critical 
point C\.j is defined over k. The normalization may be imposed at no cost if (after 
permutation of the indices) /U^i = ^ 2 ,i = M3,i = 1- 

We are now ready to detail the lifting algorithm. Out of coherent data in F p and 
a parameter N, it computes a p _w -approximation of the corresponding coherent 
data in Z p , in the form of an approximation in Z/p N . 

Algorithm 5.2 (Lift a solution p-adically) . 

Given: coherent data W t = YTjUW^f € F p [x], A e and Q l e P 1 (F P ); a 
parameter JVeN; and lifts Qi € P 1 (Z/p N ) of the points 

Return for infinitely many p: coherent data Wi = Yl^Li^i}' 3 € (Z/p N )[x] and 
A e (Z/p N ) x that reduce mod p to 

First, we assume that the data may be normalised. This amounts to requiring 
at least three of the Wij, for distinct i's, to have a linear factor. This holds for a 
positive proportion of primes p. If no such three factors exist, the algorithm aborts. 
Otherwise, we silently replace the three corresponding ft^y by {PjjjPi'j 3 in tnc 
shapes so as to create a term with /ijj = 1. 
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We write now each Wjj in the form 
(4) W id = + ]T 



w 

s=l 

N 



Hi 



for unknowns U>i,j, s € Z/p 

Recall from Q the expressions F; = W i +\q i W 1 -W 2 and J 7 = (F 3 , . . . , F k ). The 
Fj are polynomials in the variables {wij yS : ^ (1, 1), (2, 1), (3, 1)}U{A}. We lift 
the coefficients of the coherent data (Wi, A) to Z, to obtain an initial parameter a — 
(roj 2 x , . . . , w® nk „ , A ). Since the original data is coherent, we have J 7 (a ) = 
(mod p). For almost all p, the Jacobian DJ 7 is invertible by Corollary 3.3 if DT is 
not invertible at a°, then we abort the algorithm. Otherwise, we apply repeatedly 
Hensel's Lemma 5.1 to obtain a solution a to J 7 (a) = (mod p N ). 



Finally, we reconstruct the polynomials Wij out of their coefficients (which are 
just coordinates of a). 



Proof of validity. The invertibility of the Jacobian was expressed in Corollary [3T3] in 
terms of the variables Cjj. This does not make any difference: here we express them 
in terms of the Wij tS , which are elementary symmetric functions of the Cjj. □ 

Example 5.3. Consider the shapes a\ — a.2 = 03 = (4 1 ,3 1 ,2 3 ). Our example 

Wi = (x- 5) 3 (x 3 + 3x 2 + 2x + 3) 2 
W 2 = x A (x + 3) 3 (x 3 - 3a; - 5) 2 
W 3 = W 2 + iWi = {x- l) 4 (x - 3) 3 (x 3 -2x- 3) 2 

gives a vector of coefficients 

a° = (^1,2,1,^1,3,1,^1,3,2,^1,3,3, 
^2,2,1, ^2,3,1, ^2,3,2, t«2,3,3, 
^3,2,1,^3,3,1,^3,3,2,^3,3,3^) 

= (-5,3,2,3, 3,-3,0,-5, -3,-2,0,-3, -4) 

with F(w ) = (mod 11). The lift is 

a 1 = (50,-41,13,25, -19,-33,19,-60, -47,11,-46,-58, 51) 

with F(wi) = (mod ll 2 ). We can continue this process inductively. Notice that 
the precision doubles in every step. 



6. Promoting a solution from Z p to a number field IK 

If the Hurwitz problem has a solution over Z that reduces to a given solution 
over F p , then Hensel lifting will hnd it after a finite number of steps. Unfortunately 
the solutions usually involve fractional coefficients, and are usually defined over a 
finite extension IK of Q. Our first goal will therefore be to determine this extension. 

Consider a degree-e extension IK of the rationals, and a € IK. Then 1, a, . . . , a e are 
linearly dependent over Q, and therefore also over Z, i.e there exists a polynomial 

P=Po+Pit + ---+p e t e 

with all pi € Z and P{a) — 0. Let now p € N be a prime such that P splits over 
Z p ; so that we may view Q(a) as a subfield of Q p . Assume also that a is invertible 
modulo p, so that we may consider a as an element of Z p . 
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Consider now a = a (mod p N ); then we have the equation 

P(a) = Pe+lP N 

which is linear in {po, ■ ■ ■ ,p e +i}- We use the LLL algorithm [10| to find small 
integer solutions to this linear equation. The default implementation uses a simple 
heuristic to guess the correct precision N and the correct extension degree: for a 
initial precision we start with extension degree e = 1 and increase e until a solution 
is found (i.e. F(a) — mod p 2N ) or the computed shortest lattice basis vector 
norm is the same for e and e — 1. If the computed vector norm did not change, we 
increase the p-adic precision. If we have a-priori knowledge about the minimum or 
maximum expected extension degree, then it can be passed to the algorithm, which 
is more likely to find quickly a solution. 

The following algorithm is described as a process that, receiving as input an 
infinite feed of ever-more-precise approximations of a p-adic number that is known 
to be algebraic, produces an infinite stream of ever- more- likely minimal polynomials 
of that p-adic number. 

Algorithm 6.1 (Convert a p-adic number to an algebraic number). 
Given: approximations, to arbitrary precision, of an algebraic number a G Z p 
Return: polynomials P(t) £ Z[t] whose likelihood converges to 1 of being the mini- 
mal polynomial of a, as the precision of a improves. 

Assume that, for each N € N, the algorithm may receive an approximation ajsr, to 
N base-p digits, of a. The element a at is represented as an integer in {0, . . . ,p N — 1}. 

Start with d = 1 and N = 1. Then, repeat the following. Consider the lattice in 
R d+1 generated by the columns of the matrix 



(p N 


-a N 


-a 2 N .. 


• -<\ 





1 














1 





u 








1 / 



Using the LLL algorithm, find a vector (Pq, Pi, ... , Pd) in the lattice, of small norm 
6fq,d- Form the polynomial P(t) = J2i=o 

If P(oat) = (modp 2Ar ), or if d > 1 and 6n,<i = ^N,d-i: t nen output P as a 
candidate polynomial. Repeat then, after having incremented d in the first case 
holds, and doubled N otherwise. 

Proof of validity. The algorithm repeatedly increases N and d. Note that the poly- 
nomials returned may have degree < d, so increasing d is harmless, and the precision 
is increased as soon as increase in maximal degree does not improve the solution. 
Let (Po, . . . , Pd) be a short lattice vector. Then this vector is 

M ■ t {P(a N )/p N ,P 1 ,... 1 P d ) 

and in particular P(oat) = (mod p N ). On the other hand, the coefficients Pi are 
small, so P is likely to be the minimal polynomial of a. □ 

Example 6.2. Considering our example ot\ — a.2 = a-3 = (4 1 , 3 1 , 2 3 ) we found 

w l,2,l = 1400834756308742009361916361765119584358776523123371526525883115012 £ Z/ll 



16 L. BARTHOLDI, X. BUFF, C. BOTHMER, AND J. KROKER 

and Pi,2,i(wi,2,i) = (mod ll 2 * 3 ) for 

Pi,2,i(t) = 39* 6 + 117t 5 + 195i 4 + 195i 3 + 141t 2 + 63t + 16. 

Higher precision values of Wx 2 1 are also zeroes of the same polynomial -Pi, 2,1 ■ We 
take this as a hint that Pi. 2,1 is indeed the minimal polynomial of the coordinate 
w\ } 2,i in the lift of our finite field solution. 

Having found the minimal polynomials Pij. s G Z[i] for all coordinates uiij. s of 
our solution vector, we determine the field IK on which they are all defined, as the 
compositum of all field extensions defined by the Pij, s - If these field extensions 
were independent, then we should just consider all zeroes in C of the Pij, s an d 
return the corresponding rational functions. 

However, in general, the field extensions will be highly dependent. To simplify 
notation, let us assume that all Pij, s are of degree e, and that IK itself is a degree- 
e extension. Then there are e possible values for each coordinate. At this stage 
we do not know how to combine these single coordinate solutions to a solution 
vector (there are e d ( fe_2 ' possible combinations). To solve this problem we use the 
following method. 

To illustrate our method, consider a, b G IK and let P a , p,, P a +b G be minimal 
polynomials of a, b, a + b respectively. Assume furthermore that P a , p, and P a +b ai" c 
of degree e, and let a\, . . . a e , b\, . . . , b e and c\ , . . . , c e respectively be approximations 
over C of the zeroes of P a , Pb and P a +b- Consider the e x e "root compatibility 
matrix" M = (M itj ) defined by 




1 if there exists k with a.; + bj 1=3 c^., 
otherwise. 



If M is a permutation matrix, then it describes which root bj should be paired with 
a-i, namely it is characterised by Mi j = 1. In this manner, all other coordinates 
are chosen, dependent on the first choice of a root of Pi 2,1 • 

Example 6.3. Tentative minimal polynomials of 1^1,3,1 and a = wi,2,i +1^1,3,1 are 

p 1 3 1 = 28431i 6 + 255879i 5 + 982449i 4 + 2056509* 3 + 2465721* 2 + 1597239t + 439138 
P a = 28431t 6 + 341172* 5 + 1844856* 4 + 5660928t 3 + 10384524t 2 + 10807344* + 5068144 



We obtain the following approximate zeroes over C using Brent's method, im- 
plemented in PARI [2]; we preserve the ordering in which the roots were returned. 



Wl,2,l 

-.150 + .807i 
-.150 -.807* 
-.5 + A40i 
-.5 - A40i 
-.850 + .807* 
-.850 - .807i 



^1,3,1 
-1.5 - 1.02i 
-1.5 + 1.02* 
-2.012 - .377* 
-2.012 + .377* 
-.988 - .377* 
-.988 + .377* 



•U+2,1 + w l,3,l 

-2.161 - 1.184* 
-2.162 + 1.184* 
-2 - 1.462* 
-2 + 1.462* 
-1.839 - 1.184* 
-1.839 + 1.184* 



ALGORITHMIC CONSTRUCTION OF HURWITZ MAPS 



17 



Now consider the compatibility matrix M . It is 



/0 1 

10 

10 

1 


\0 1 







1 



0/ 



i.e. Mij is 1 precisely when (u>i,2.i)i + (wi,3.i)j approximates one of the values 
in the last column of the table above. We note that M is a permutation matrix. 
Applying the permutation M to the list of roots (wi j3j i), of -Pi. 3,1 leads to a valid 
coordinate pairing: now {Mw\^^) i corresponds to the (wi,2,i)r 

Example 6.4. There are situations in which the matrix M is not a permutation 
matrix but nevertheless contains useful information. Consider the finite algebraic 
set 



In this case the minimal polynomials of the coordinates are x 2 + 1 and y 2 — 2 and 
the root compatibility matrix is M — ( \ \ ) . All 4 possible pairings lead to correct 
solutions. 

There are also situations in which the matrix M differs from the permutation 
matrix giving the correct identification of coordinates, even with exact arithmetic. 
For example, if £ denotes a fifth root of unity, and 



then the coordinates have respectively x 4 + x 3 + x 2 + x + 1 and y 4 + y 3 + y 2 + y + 1 
as their minimal polynomial, and the root compatibility matrix is insufficient to 
recover S. Indeed the sum of the coordinates z = x + y, which has minimal 
polynomial z 4 + 2z 3 + Az 2 +32 + 1, docs not distinguish solutions in S from the 
non-solutions 



Note that S and 5" are distinguished by the equation x 2 — y which holds in S but 
not in S' . A linear form ax + by with O^a^i^O would also do. 

There are also situations with more than two variables in which all compatibil- 
ity matrices between two variables lead to possible pairings, but their combined 
information is not enough to identify the correct solutions. For example, consider 



All pairings between two coordinates are allowed, but there are 4 solutions in total, 
and not 8. 

Above we have used the linear forms Xi + Xj to determine the compatibility. Our 
algorithm uses random linear forms to avoid these problems, or at least make them 
less probable. 

We are now ready to explain the algorithm computing the solutions in number 
fields that reduce modulo p N to given approximate solutions in Z p . That problem 
is in fact an instance of the following, more general problem. 



S = { (*, V2) , ( - i, V2) , (*, -V2) , ( - i, -V2) } . 



5' = {(e,C 3 ),(C 2 ,0,(e 3 ,C 4 ),(C 4 ,e 2 )} 
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The following algorithm is described as a process that, receiving as input an 
algebraic system over Z and an infinite feed of ever-more-precise p-adic approxi- 
mations of a solution, produces a stream of algebraic solutions (in the form of a 
minimal polynomials over Z and a complex numbers singling out roots of the min- 
imal polynomials). The stream eventually exhausts all solutions conjugate to the 
p-adic solution. 

Algorithm 6.5 (Solve O-dimensional algebraic systems). 

Given: a polynomial system of equations T = (Fi, . . . , F m ) in variables x%, . . . , x m , 
having a finite number of solutions; and approximations, to arbitrary precision, of 
a solution (x~i, . . . , x^) in Z p 

Return: a number field K c C, and exact solutions (x\, . . . , x^) in IK, for i = 1, . . . , s; 
each element of IK is given by its minimal polynomial and an approximation in C of 
a particular root. 

We construct the solutions (x\, . . . , x l m ) iteratively, entry by entry, by construct- 
ing partial tables {{x\, . . . ,xV) : i = 1, . . . , s}. We start by an empty table, with 
s = 1 and k = 0, and let e be a small number. 



Then, for each k = l,...,m, we do the following. Using Algorithm 6.1 we 
compute a likely minimal polynomial of x~k, say of degree e. We compute, to 
precision better than e, approximate roots r±, . . . , r e € C of P^. Set 

Si = min{|r; — rj \ : 1 < i =/= j < e}. 

If < e, we halve e and restart all over. 

We next choose randomly a linear form L = liXi + • • ■ + IkXk with ^ 0. Again 



using Algorithm 6.1 we compute a minimal polynomial Pl for L(x~i, . . . , x~k). If the 
degree of Pl is not divisible by s, we choose a different linear form L and repeat 
the above. Otherwise, we let Si be the minimal distance between roots of Pl. If 
5i < e, we halve e and restart all over. 

We then compute the s x e matrix M — (Mij) with 



ii\P L (L(xi,...,xi_i, rj ))\ < !mm{<M2}||£|| 
otherwise. 



Let the degree of Pl be st. If M contains t ones per row and one one per column, 
then we replace s by st and replace each row (x\, . . . , x l k _ 1 ) is the partial table by 
t rows (x\, . . . , x\ rj) for all j such that Mij = 1. Otherwise, we repeat the 
step with a different linear form or, if that failed more than ten times in a row, we 
simply skip the iteration. 

When the iteration finished with k = m, we have obtained s candidate solutions, 
which we check algebraically by evaluating T on them. We output all those that 
are certifiably valid solutions, and restart the algorithm with better approximations 
of the xi, . . . , x^,. 

Proof of validity. First, all the solutions returned are valid, since they were checked 
(using exact algebra) by evaluating J- on them. 

Let now (xi, . . . , x rn ) be a solution that is conjugate to (a?T, . . . , x^). In particu- 
lar, the minimal polynomials of the Xi and x\ are the same, so they will eventually be 
found by Algorithm |6.1| Similarly, for every linear form L with integer coefficients, 
L(xi, . . . , x m ) and L(x~i, . . . , x^) also have the same minimal polynomial, so it will 
also be eventually found by Algorithm |6.1| □ 
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We apply this algorithm to the same polynomial equations (|3| . The variables Xi 
are a relabeling of the u>i,j, s from Q. 

7. Computing the monodromy 

In this section, we detail the second part of the algorithm sketched in f|2] 
We are given an approximation of a degree-d rational map / G C(z), as well as 
an approximation of the critical values Q C C, and the local degrees Q(g,i, • ■ ■ , otq t t 
above each critical value. We are asked to compute the monodromy of the covering 
induced by /. 

The first step is to compute a triangulation =2 of P : (C) by arcs of circle, and 
containing Q among its vertices. A particularly efficient triangulation is the Delau- 
nay triangulation. This is a decomposition of P 1 (C) into triangles, such that, for 
any two triangles with a common edge, the sum of their opposite angles is > tt. 
Such a triangulation always exists; is essentially unique; and may be computed e.g. 
using [l4] . 

For performance reasons, we refine the triangulation by adding vertices to it: 
whenever we encounter a triangle whose ratio "circumradius / shortest side" is 
larger than 1000, we add the circumcenter to the triangulation. This process con- 
verges, and gives a reasonably good triangulation in that its triangles are not too 
acute; see [15] , 

The dual decomposition J2 1 - of the sphere is the associated Voronoi diagram. 
It has one vertex, called a dual vertex, per Delaunay triangle and one edge, called 
dual edge, across every Delaunay edge. Each of its edges e is parametrised as the 
prcimage, under a Mobius transformation of the arc [0, 1]. We denote by W 
the vertex set of J2 1 - , and choose a basepoint * € W . For each w £ W, we number 
arbitrarily the elements of the fibre f~ 1 (w) as {wi, . . . , Wd}- Because W is far from 
Q, there are d preimages of each w G W , and their computation is numerically 
stable. 

There are now two strategies, which have both been tested and implemented. 
The first one is a bit simpler, but the second one performs better in practice. Both 
associate a permutation ^ e with each edge e € , in such a way that the the /-lift 
of e that starts at w'j ends at w" qn. Both are explained below; assuming them, we 
finish the description of the algorithm. 

For each critical value q 6 Q, let e(l), . . . , e(n) be a path in J2 1 - that starts and 
ends in *, and surrounds once counterclockwise the point q and no other vertex of 
Q. These paths form a basis for the fundamental group cri(P 1 (C) \ Q, *). Compute 
the permutation a q = <; s m • ■ • S £ (n)- Then the monodromy representation of / is 
given by the family (a q ) qeQ . 

7.1. Connect-the-dots. For each dual edge e G £i , going from w' to w" , we do 
the following. Knowing the spherical distance from w' to w" and using coarse esti- 
mates on |/'(z)|, we have an upper bound on the length of each of the d preimages 
of e. We attempt to match each w'j with a , w",^ for some permutation ? e 6d, 
by matching each Wj to the closest w"uy If more than one match is possible, we 
subdivide the path e; at some time we will have selected points e\, e M on e for some 
< A < jj, < 1, and we may need to insert a midpoint £(a+ m )/2 between them. We 
then compute its d preimages and attempt to match them to the e 3 x and e*. After 
sufficient subdivision, we have obtained d lifts of e, which are all well separated 
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from each other, and we know that the lift of e that starts at Wj ends at w"tj\ for 
some permutation q = <; £ . 

7.2. Using two triangulations. We also compute the Delaunay triangulation 
on f^ 1 (Q), and parametrise its edges e via Mobius transformations v e such that 
e = i/ f - 1 ([0,l]). 

It is straightforward to lift £&^~ through /: its edges are all the curves defined by 
equations (fj, e o f)(z) € [0, 1]. We consider as before an edge e G =2^ , going from 
w' to w". For each i = 1, . . . , d, we consider the starting point u^, and wish to find 
j G {1, . . . , d} such that the lift of e starting at w[ ends at w'J. 

We first determine in which triangle T of ^ the lift w\ lies. Then we compute 
whether the lift e of e starting at w[ leaves T. If this happens, then it must cross an 
edge e of T, namely, we have /i e o/oj/~ 1 ([0, 1]) € [0, 1]. This entails, firstly, that the 
imaginary part of fi e o / o /j,^ 1 vanishes, and secondly that its real part belongs to 
[0, 1]. Both arc polynomial conditions imposed on real- valued polynomials, and are 
efficiently computable numerically. We also keep track of the point of intersection 
of i and e. 

In that case, we move to the neighbouring triangle T" of T along edge e, and 
continue. When we do not detect more intersections with edges of , we know in 
which triangle of the vertex w'J lies. 

It may happen that two or more vertices w'J belong to the same triangle T that 
we have found in the previous paragraph. In that case, we let w denote the last 
point on i that was computed — possibly w[; it also belongs to T. We consider 
in turn all candidates w'J , and compute the straight path 5j from w to w'J and its 
image f(Sj). If there exists a unique j such that f(Sj) lies in the two triangles 
of £2 to which e belongs, then we have found the desired j. If there are no such 
j, then we interpolate. Writing w = e(to), we consider t € (to, 1), and those lifts 
w e f^ 1 (e(t)) that belong to T; we then consider the paths 5j from w to w'J as 
before, and continue with increasing t. 

Additional care must be taken for points that lie on edges, and not inside trian- 
gles, of the triangulations; and of crossings of edges with multiplicities. However, 
because of the necessary crudeness of norm estimates on |/'(z)| in the first method, 
this second method is preferable. 

8. AN APPLICATION TO DYNAMICAL SYSTEMS 

The post-critical set of a branched self-covering / : S — > S with branch value set 
Qf is 

p f ■■= U / on (Q/)- 

n>0 

We are interested in the case where Pf is finite and we consider / up to isotopy 
rel Pf, namely, / ~ g if there exists a path of branched self-coverings from / to g 
whose post-critical set moves smoothly. We say that / is equivalent to a rational 
map F if there are homcomorphisms <j> : S ^ P 1 (C) and ip : S — > P 1 (C) such that 
F o (p = tjj o / and (j) is isotopic to ip rel Pf . 

On the one hand, many examples of branched self-coverings can be constructed 
combinatorially, via triangulations; for these, it is natural to consider the maps up 
to isotopy. On the other hand, a fundamental theorem by Thurston claims 
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Theorem 8.1 (Thurston; see [6]). Let f : S — > S be branched self-covering with 
jfPf > 3. Then f is equivalent to a rational map if and only if f admits no 
"Thurston obstruction" , namely, if and only if, for every collection C of isotopy 
classes of non-peripheral curves on S\Pf, the QC-endomorphism 

C 3 c h-> V - — , ■ d e QC 

deg(/ : d -> c) 

/las spectral radius < 1. 

Furthermore, in that case, the rational map is unique up to conjugation by a 
Mobius transformation. 

8.1. Cui's Problem. Cui Guizhen suggested in 2010 that a "Sierpinski map", 
namely a rational map whose Julia set is a Sierpinski carpet, should have an in- 
variant non-peripheral curve with unbounded number of preimages. He then found 
a counterexample, given combinatorially as follows: 




In that case, Pf = Q := {Qi, Q2, Q3} with Qi := 00, Q2 := and Q 3 := 1. From 
this picture, it is easy to compute the monodromy action about Q: the permutations 
are those given in |2]), namely 

a 1= (1, 7, 11, 2)(3, 8)(4, 5)(6, 10) (9, 12, 13), 

<r a =(l,3,12, 4) (5, 9) (6, 7) (10, 13, 11) (2, 8), 

(73= (1, 5, 13, 6)(7, 10)(2, 3)(8, 11, 12)(4, 9). 

Because ffQ = 3, all curves on P 1 (C) \ Q are peripheral. By Thurston's rigidity 
theorem, there is a unique map with monodromy 02, 03) and fixing 00, 0,1 
with local degree 2, so it is the map computed by our algorithm (perhaps after 
precomposition with a Mobius transformation so that 00, 0, 1 are critical of order 
2). 

Our algorithm searched in fact for a map with 00, 0, 1 of order 4. This is an 
improvement to searching immediately for the correct map, because there are three 
points of order 2 above each of 00, 0, 1, and they may lie in a strict field extension. 
Needless to say, it suffices, once such a map is found, to precompose it with a finite 
number of possible Mobius transformations. 

Our algorithm found a solution ( mod 11) of the defining equations for a map; 
then lifted them ( mod ll 2 ) and finally obtained six Galois conjugate solutions. 
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The correct one (with correct choice of point of order 2 above oo,0, 1) was then 
found. The original map is of the form 

fiz/w) - A — - hiW ^^ z ~ b z w ) 3 ( z ~ h 2,iw) 2 {z - b 2 ,2w) 2 {z - b 2 , 3 w) 2 , 

here are the preimages <ij, 6j, Cj of oo, 0, 1 respectively: 

04 = 00 (meaning the term z — a^w should be replaced by 1), 
a 3 « 0.500000000000000000000000000000 - 0.439846359796987134487167714627*', 
a 2 ,i ~ 1.61268567872451072013417667720 - 0.490182463946729812334860743821*', 
a 2 , 2 ~ 0.500000000000000000000000000000 - 0.0415300696430258467988035191529*, 
a 2 , 3 ~ -0.612685678724510720134176677204 - 0.490182463946729812334860743821*, 
b 4 = 0, 

b 3 w 1.12748515145901194873474709466- 0.991840479188802206853242764751*', 

62.1 ~ 1.98629656633071582984701575517- 0.164982069462835473582606346591*', 

62.2 ~ 0.567640411622375679553529964298 - 0.172536644477962176299255320022*', 
b 2 , 3 « -0.995164705141609432502666361446 - 0.796186860797306011242450678339*', 

c 4 = 1, 

c 3 w -0.127485151459011948734747094655 - 0.991840479188802206853242764751i, 
c 2 ,i w 0.432359588377624320446470035702 - 0.172536644477962176299255320022i, 
c 2 ,2 « -0.986296566330715829847015755165 - 0.164982069462835473582606346591i, 
c 2 ,3 « 1.99516470514160943250266636145 - 0.796186860797306011242450678339*, 

A w 0.130027094895701439414281708196*. 

The required Mobius transformation /* maps (02,2:^2,3122,2) to (00, 0,1) respec- 
tively, so 

, , Z- 6 2 ,3 C 2 , 2 - 0,2,2 
K z ) = 7 

z - a 2 , 2 c 2 , 2 - 2 ,3 
and the required solution is / o /*. Its Julia set is displayed in Figure [T] 
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